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Abstract 

We study the collapse of a circular loop of cosmic string. The gravitational 
field of the string is treated using the weak field approximation. The gravi- 
tational radiation from the loop is evaluated numerically. The metric of the 

loop near the point of collapse is found analytically. 
02.60.Cb, 04.30.+X, 98.80.Cq 
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I. INTRODUCTION 



In the cosmic string scenario |0 a network of strings forms in the early universe and 
evolves. This network contains loops of various sizes and shapes. Typically a loop oscillates 
between some minimum and maximum size. The oscillation produces gravitational radiation 
and the loop gradually shrinks and disappears. If, however as the loop oscillates its minimum 
size becomes smaller than its Schwarzschild radius, one would expect the loop to collapse 
to form a black hole with a concomitant emission of gravitational radiation. In particular, 
a circular loop, whose minimum radius is zero, would be expected to form a black hole. 
Now in the case of a circular loop of cosmic string there are two effects which could be of 
significant influence in the dynmaics of the gravitational field that would be produced. We 
take the point of view that up to times at which the loop radius is a few orders of magnitude 
larger than the Schwarzschild radius, the weak field approximation for the metric should be 
a reliable guide to the character of the fields produced. As the loop radius monotonically 
decreases, its speed quickly approaches the speed of light. Thus one would expect that 
regions inside the loop would not have been affected by the data on or outside the loop for 
such late times. However, there is another aspect of the dynamics in operation here. The 
loop, while quickly approaching light speed, is also accelerating. One would thus expect the 
world-sheet to exhibit a build-up of larger and larger fields in the vicinity of the loop. In 
the collapsing loop problem, it would seem that both effects are in operation and neglecting 
one in favor of the other would be a dubious move. Hawking has argued that a circular 
loop does form a black hole and that no more than 29% of the mass of the loop is lost in 
gravitational radiation. His argument is based on the approximation that the loop moves 
at the speed of light and that the metric of the interior of the loop is fiat. 

In this paper we treat the collapse of a circular loop in the approximation of weak field 
gravity. In section II we set up the equations for the metric of the loop. Section HI contains 
a numerical evaluation of the gravitational radiation produced by the loop, including the 
power radiated per solid angle and the total power radiated up to late times. In section 
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IV we find an analytic expression for tlie tidal force due to the loop in the limit as the 
loop gets small. We show that the interior of the loop is fiat. We also find the late-time 
weak field metric compatible with the tidal force for the region exterior to but nearby a 
critical cylinder where the fields exhibit a large spatial variation. Section V contains our 
conclusions. Throughout this paper we use units in which c = G = 1. 



II. METRIC OF THE LOOP 

In Minkowski spacetime the motion of a loop of string is specified by x'^ = f^^{(T,T), 
that is by giving the spacetime coordinates of the string as functions of two worldsheet 
parameters a and r. With the usual choice [§] of the parameters a and r the string equation 
of motion can be shown to be — f"^ = 0. Here a dot denotes derivative with respect to 
T and a prime is derivative with respect to a. The solution of this equation for a circular 
loop is /° = r and 

-* T / a a \ 

/ = — ro sin — [ x cos — + y sin — . (1) 

'^o V ro ro/ 

This is a loop that begins at rest with radius rg at time — 7rro/2 and collapses down to a 
radius of zero at time zero. 

The stress-energy of the string is 

T'^^ = /i Z^'" da {rr - rr) ^^^^ - /(a, t)) (2) 

where /i is the mass per unit length of the string. 

In the weak field approximation the metric is rj^j^y + h^j_y where rj^y is the fiat metric and 
/i^i/ is small. Einstein's equation is then 

where S^j^y = T^y — | Trjfj_y. In obtaining equation (3) we have also used the harmonic 



coordinate condition d^h^y = where hny = h^y — Ihrj^y. Suppose initial data for h^y i 
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given on the surface t = to- Then the solution of equation (3) is h^y = H^y + I^y where H^y 



is the part due to the stress-energy and 1^,^ is the part due to the initial data. ||^ (see Fig. 

First consider H^y. Let be the point at which we evaluate H^^. Let be a point 
on the worldsheet of the string and on the past light cone of x^. Define and W by 
k'^ = x^ — and W = — /c^/^. Then if^j, is given by 

^M. = 4/i£"rfa^ . (4) 

Here F^,^ is given by F^^ = f,i fu - f/i' fJ - rj^y and the integral is over all points 
in the past light cone of x^ and to the future of the initial data surface. Equation (4) 
can be derived using the Green's function for the wave equation and the form of the string 
stress-energy tensor. Define the angle A by A = (a/ro) — 0. Then equation (4) becomes 

H,. = 4^MroJ^^ dX^ . (5) 

Since is in the past light cone of x'^ we have k^k^ = Q. It then follows from equation 
(1) that 



0= r -r + rosin^ + 2i(;r + 2r sin 6^ ro sin — cos A . (6) 

^ ^ V ro / ro 

Similarly from equation (1) it follows that W is given by 

T T T 

W = t — T + r COS — sin 6 cos A + ro sin — cos — . (7) 

ro ro ro 

Equation (6) gives A as a function of r or gives r as an implicit function of A. This result 
will be used in the next section to numerically integrate equation (5) and in section IV for 
an analytical treatment in the limit as the loop gets small. 

Now consider J^j,^. Let S be the intersection of the past light cone of x^ with the initial 
data surface. Then 5* is a sphere of radius R = t — tQ. Define P^,^ (r) and V^i, (r) by 
P/iu (r) = h^^ {to, r) and Vf,u (r) = h^^ {to, r). Then we have 

V = ^ ^ [R''P^.u + n-VP^, + V^,] R^dn . (8) 



Here the natural volume element on S is used and n is the outward pointing normal to 5*. 

To find 1^1^ we must first choose the initial data P^^, and V^j,^. We choose the time to to 
be — 7rro/2, when the loop is at rest. Since to is a moment of time symmetry for the loop we 
choose Vfj^u = 0. The initial data must be chosen consistent with the harmonic coordinate 
condition d^h^y = 0. We choose P^j^ to be the solution of V^P^^jy = — IGvrS'^y. It then follows 
that I^i, is given by 

V(r,t) = P,.(f) - 4 / d'r' . (9) 

Jv \r — r'\ 

Here V is the interior of the sphere S. 

III. GRAVITATIONAL RADIATION 

We now calculate the gravitaional radiation produced by the loop. Define the retarded 
time u hj u = t ~ r. Then the radiation zone is given by the limit as r — oo with u 
remaining finite. All calculations in this section will be done in the radiation zone. The 
power radiated per unit solid angle in gravitational radiation is given by 

dP _ ( dh"'^ dKp 1 
'dn~mi I du ~du 2 

In the radiation zone the harmonic coordinate condition becomes dh'^^/du = 0. Since the 
loop is not rotating we have h^[adf3](f) = 0. It follows from these two conditions that 

dP _ J_ ( dh^^/duV 

dn ~ IQtt \ rsin^e J ' ^ ^ 

We therefore need to find H^^ and I^^j, in order to find the power radiated. 

Define ip and w hj ip = t/tq and w = u/tq. The radiation zone limit of equations (6) 
and (7) yields 

w = ^p + smip sin 9 cos A , (12) 
W = r {1 + cos 'i/' sin 6^ cos A) . (13) 



dh 
du 



(10) 
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Prom equation(l) we find that F^^ — r^sin ^^sin^A. Now define B hy B'^ = I + 
cos'0 sin^ cos A. Then we have W — rB~^ and thus 

^'"'^ 4/iro dASsin^A . (14) 



rsin^^ 



The range of integration for A is restricted by the fact that > — 7r/2. Thus A± are 
given by 

TT 

— + w + sin^ cosA± = . (15) 

If (7r/2) + w > sin 6* then the range of A is from to 27r. 

We now find /<^<^. Taking the radiation zone hmit of equation (9) and using equation (2) 
we find 

14,4,= P^ir) - Ajirorsin'^O / dA sin^A . (16) 

v A— 



It then follows that 



rsin 9 rsin 9 



4//ro f^^ dX {B - 1) sin^A . (17) 

J A — 



We can now compute the power radiated. Taking the derivative of equation (17) with 

respect to u we find 

dh4>4>/du _ f>^+ dB .^2 



4u / dX — sin^A . (18) 
rsin'' 6* Jx- dw 

(Note: though X± depend on u, the quantity B — 1 vanishes at A = A±. Thus there is no 
dependence on dX±/du in equation (18) ). 

Equation (12) gives A as a function of ip; however, this equation gives ip only as an 
implicit function of A. We therefore change the variable of integration from A to ijj. Define 
the integral K by 

dX 



K 



Jib- 



dip 

Then since the range of A covers the range of -0 twice we have 
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sin^A^ . (19) 
dw 




(20) 



Expressing all quantities in equation (19) as functions of ip we find 



K 



(21) 



V'- sin 9 [sinip + {w — ip) cosip] 



The range of integration is fixed by two conditions: first ip must lie in the interval 
[— 7r/2,0]. Second the argument of the square root in equation (21) must be positive. We 
numerically evaluate the integral as follows: for each value of w and 9 we find the allowed 
range of ip. We then numerically integrate to find the integral K and use equation (20) to 
find the power radiated. 

The results are shown in Fig. |^ and ^. Fig.|^ shows the power radiated as a function 
of w and 9. Fig. | shows the total power radiated as a function of w. Note that most of 
the power is radiated near w = i.e. when the loop is small. This is because the loop 
starts out at rest and then accelerates inward. As w approaches zero the speed of the loop 
approaches the speed of light. The large speed and acceleration of the loop at late stages 
of the collapse gives rise to a large amount of radiation. Also note that most of the power 
is radiated near 9 = 11/2 i.e. near the plane of the loop. As the loop shrinks the power 
becomes more narrowly peaked around this plane. 



The loop collapses down to a point at t = 0. We now find the gravitational field at 
spacetime points near the point of collapse. That is we consider values of r and t much 
smaller than tq. For the weak field approximation to be valid the loop must be much larger 
than its Schwarzschild radius. That is we must have \t\ » ^tq. For the cosmic strings 
considered in cosmology /x ~ 10^^. Thus there is a range of t, /xro << \t\ « tq where 
the approximation of this section is valid. We implement this approximation by considering 
t and r to be of order 1 and taking the rg ^ 00 limit. In this limit we have t/tq « 1. 



IV. FIELDS NEAR THE POINT OF COLLAPSE 
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However, as we will see, r can be much larger than t or r. So we will have to keep lowest 
order nonvanishing terms in t/tq. 

In the large tq limit the tensor I^i, vanishes. It then follows that the metric perturbation 
is given by 

2m F„ 



Here m — 27rro// is the mass of the loop. 

We now examine the behavior of W. Taking the large Tq limit of equations (6) and (7) 
we find 

/ N 1 

O^(r^-f) + 2r (t + rsin^cosA) , (23) 

^ ^ 3 Tq 



2 

= i + rsin^cosA ^ . (24) 

3 



Defining U by U = t + rsinO cos A and using the last two equations to eliminate r we find 

16 (r^ - t ^f 
3 



{W + Wf {W -U)^^ ^ 2^ . (25) 



Note that the right hand side of this equation is small. Thus we have W ^ U or W ^ —3U. 

When U passes through zero W gets very small (in fact it approaches zero in the ro — > oo 
limit.) This makes diverge in the limit. However, it turns out that this divergence is 
gauge. Physical quantities remain finite. We will therefore work with the Riemann tensor 
which is gauge invariant. The Riemann tensor is given by 

-R/ui/a/3 = df3d[nhi,]a — dad[^hvT^j3 . (26) 

To find the Riemann tensor we begin by calculating d^hfj,,^. 

2m /F„„^ F, 



dzh 



Z"'(J,V 



Prom dz{k^k^j) = we find d^T = —z/W. Then from d^W = — ^^(A;^/'^) we find 

= ^ {kj' - rU) ■ (28) 
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Evaluating the quantity in parentheses in equation (28) we obtain 

_-2m r^^dX (F,, F, f rsingcosA \\ 

We also need to compute h^z- In the large tq limit F^z — T^/rg so it follows that 

2m dX 

Define the critical cylinder to be the surface where r sin ^ = —t. Then the behavior of 
the metric inside the critical cylinder is very different from the behavior outside. 



A. Inside the Critical Cylinder 

We now consider the case rsin^ < —t, that is inside the critical cylinder. Then we have 
U < and U is bounded away from zero. Therefore W is bounded away from zero. It then 
follows that in the large Tq limit we have hzz 0. Furthermore since F^^ — > in this limit it 
follows that dzhfj,!, — > 0. It then follows that Rabcz — 0. However, the spacetime off the string 
is vacuum. So the Riemann tensor is equal to the Weyl tensor and so Cabcz — 0. But it 
follows from the properties of the Weyl tensor that if Cabcz = then the Weyl tensor is zero. 
So in our case the Riemann tensor vanishes inside the critical cylinder and the spacetime in 
this region is flat. 



B. Outside the Critical Cylinder 

Now we treat the case where rsin^ > —t, that is outside the critical cylinder. Deflne 
the variable X by 

2(^2 _ ^2)3/4 

Expressing W, U and r in terms of X we have 
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2^/^ V 3 
t2 Vr2^ 



2 [X - , (33) 



(34) 



Prom these equations we can express all quantities in terms of X and Tq. 

We now compute the metric component hzz- The quantity F^z/W is a bounded function 

— 1/2 

of X multiplied by Tq . It then follows that in the large Tq limit h^z 0. 

Now we compute dzh^j,^. In the large Tq limit the only contribution to the integral comes 
from the region where W is very small. That is the region where U ~ This is the region 
where A ~ Aq where Aq is given by t + r sin 9 cos Aq = 0. In the integral for dzh^i, we change 
the variable of integration from A to X. Since A ~ Aq in the region of interest we should 
be able to evaluate all quantities at A = Aq. Define the coordinates p and s by p = rsin^ 
and s = v^p^ — t^. Then from equations (32) and (33) it follows that dU/W = dX/X 
and therefore that dX/W — dX/{sX). Our approximation that A ~ Aq is valid provided 
that the fractional change in X is much larger than the fractional change in sin A in the 
region where U ^ 0. That is we require \dX/X\ » \dsmX/ sinA|. This is equivalent to 

» \t\{r'^ — t'^)^^'^ / y/r^. This will be satisfied except for extremely small values of s. 

Now the range of A covers the range of X twice. So for functions that are odd in A the 
integral vanishes. For functions that arc even in A there is a factor of two multiplying the 
integral. We then find that dzh^i, is given by 

Now express all quantities in the integrand in terms of X and tq. Then we find that F^i, 
is of order 1, is of order Tq ^, T/{Wro) is of order 1 and t/r is of order 1/y^. So 

the F^i, and the t/2T terms in the integrand give zero contributions in the large ro limit. 
Define F^^, by FfMu = Ff^i,\ . Then since the integrand is peaked around Aq and since F^i, 
is slowly varying we have 

_ -8m zF^, AX^dX 
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In the large tq limit the limits of integration are zero and infinity. So we find 

d.K. = ^ .fr\. ■ (37) 

Evaluating Ff^^y at A = Aq we find 

Ff,. = d,{t/p) d,{t/p) + d,<P 9,0 . (38) 

We now compute the Riemann tensor. Since h^z = it follows that 

RfiuTz = ^ {d^dzKr - d^dz^h^r) ■ (39) 

So R^uTz can be computed using the expression in equation (37). However, in vacuum 
the Riemann tensor is equal to the Weyl tensor and is determined by R^urz- Let and 
z'^ be the unit vectors in the and z directions respectively. Define the tensor by 
E^i, = C^^,y^- Because the axial Killing vector d/dcp is hypersurface orthogonal, the Weyl 
tensor is determined by E^^. We have 

Cf,ura = -40[^^,][^0^] - e^^'^Pera'''^E^n4>p(j)q ■ (40) 

The tensor Efj,i, is determined by Cf^i^rz by 

Ef,^ = (0/,0, - Vi^u) - C^zu, - 2i(^C^)^^^ . (41) 

Finally some straightforward but tedious algebra yields the following expression for E^^. 
4m s 



"^/^^ ^ 9 / 9 , 9n2 



(s^ - z^) [d^zd^z - df^sdys) - 4zs d(^zdy)S . (42) 



This expression is actually singular as s — > 0. (that is as p — —t). This is because is 
singular as s — > 0. {dpS = p/s). So the Cartesian components of the Riemann tensor diverge 
like 1 / y/p + t as p ^ — t. Fig. |^ shows the Riemann tensor component Rp^p;j, as a function 
of p and z. The component Rp^pp^j, is measured in units of Am/nt^ and the coordinates p and 
z are measured in units of |t|. 

It turns out that the quadratic curvature scalars are not singular (except at 2; = s = 0). 
We have 
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8m ' ^ 



(43) 



Note that things are becoming singular just where our approximation breaks down (i.e. 
for extremely small p -\-t). However, for moderately small p -\-t our approximation is still 
good. So we find that the tidal force shows rapid variation and large values near the critical 
value of p. 

The metric perturbation /i^,^ diverges at late times. However, except for the critical 
cylinder this bad behavior is pure gauge. This means that there is a well behaved metric 
perturbation that gives the same Riemann tensor as h^^. Define the perturbation h^^ by 

4m - 



h 



— F,^s-Hn(s' + z') . (44) 



Then it follows from equation (37) that dzh^y = dzh^i,. It then follows from equation (38) 
that h — and d'^hfj^v = so h^i, is in harmonic gauge. Finally some straightforward but 
tedious algebra using equation (38) shows that d^d^h-rrr — 0. Therefore h^i, is a solution 
of the linearized Einstein equations. However, since dzh^i^ = dzhf^y it then follows that the 
Riemann tensor of hfj_i, is the same as that of h^^. Therefore h^i, can be regarded as the 
metric perturbation of the string loop at late times exterior to the critical cylinder. 



V. CONCLUSION 

The calculations of the last two sections give a picture of the fields produced by the 
collapse of a circular loop of cosmic string. When the loop begins to collapse it is moving 
slowly and its gravitational field is quite mild. Very little power is radiated in gravitational 
radiation in this early stage of collapse. As the loop gets smaller it moves faster and more 
power is radiated. Most of this power is radiated in directions near the plane of the loop. 

When the loop gets very small its speed approaches the velocity of light. Near the loop 
the tidal forces become large. These tidal forces have a strong spatial dependence. We have 
shown that inside the critical cyhnder the tidal force vanishes so the metric is flat. This 
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flatness of tlie metric was assumed by Hawking in liis argument tliat tlie loop must form 
a black hole. Thus we have supplied a missing step in Hawking's argument. 

Outside the critical cylinder there is tidal force and the tidal force becomes large as one 
approaches the critical cylinder (strictly speaking calculation gives an infinite tidal force on 
the critical cylinder. However our approximation breaks down when we get too close to the 
critical cylinder. Thus all we can say is that near the critical cylinder the tidal force has 
large values and large spatial variation.) 

When the size of the loop gets to ~ 1/ ^ oi its original size it has approached its 
Schwarzschild radius. In this regime the weak field approximation, which we have used 
throughout this paper, breaks down. Thus our results do not tell us what happens when 
the loop collapses to form a black hole. To treat this strong field regime one would need to 
use general relativity rather than the weak field approximation. The equations then become 
much more complicated, and there is little hope of an analytic treatment. A numerical treat- 
ment may be feasible. The results of section IV could be used as initial data for a general 
relativistic numerical evolution of the string and its gravitational field. Note, however, that 
the fields have large spatial variation. Thus one would need a computer code that could 
handle large spatial changes in order to evolve the loop until it collapses completely to form 
a black hole. 
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FIGURES 

FIG. 1. Spacetime diagram showing the intersection of the past hghtcone of the field point and 

the world-sheet of the collapsing string loop. 



FIG. 2. Plot of power radiated per unit solid angle as a function of w and 9. 

FIG. 3. Plot of the total power emitted as a function of w. 

FIG. 4. Surface plot of Epp as a function of z and p. Epp is measured in units of Am/Trt^ while 
p and z are measured in units of \t\. 
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